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Q Abstract 

^sO The fluid dynamics video that is presented here outlines recent advances in the simulation 

of multiphase cellular blood flow through the direct numerical simulations of deformable red 

1 blood cells (RBCs) demonstrated through several numerical experiments. Videos show particle 
deformation, shear stress on the particle surface, and the formation of particle clusters in both 
Hagen-Poiseuille and shear flow. 
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Introduction 

This video includes demonstrations of dilute and dense suspensions of RBCs in wall-bounded shear 
and pressure-gradient driven pipe flow. The video, available in hi-resolution and web-friendly 
formats, can be found at http://ecommons.library.cornell.edu/handle/1813/14097 

The animations presented in the movie submission are results from the direct numerical sim- 
ulation of deformable suspensions using a hybrid lattice-Boltzmann/finite element method that is 
coupled with a particle-particle and particle-boundary contact model [1] . The code has also been 
optimized to run on as many as 65,536 cores of the IBM Blue Gene/P architecture [2]. In order to 
simulate high volume fraction (hematocrit) simulations of deformable particles in wall-bounded or 
cylindrical domains, a seeding method is used. In this procedure, the particles are initialized at 30% 
of their final size and grown to their full size rigidly, allowing particle-particle and particle-wall 
interactions using the lubrication and contact model previously referenced. Each RBC is composed 
of 504 linear triangular shell elements (254 nodes) resulting in a 762x762 finite element matrix for 
each particle updated with Newmark's time integration method. Physically, RBCs are biconcave in 
shape with a maximum diameter of 8 urn enclosed by an elastic membrane with an effective elastic 
J> shear modulus of 5.7 x 10 _3 dyn cm -1 [3] and a bending stiffness of 2.2 x 10~ 12 dyn cm [4], which 

surrounds a liquid hemoglobin with a viscosity of 6 cP. The RBCs are suspended in blood plasma 
with a viscosity of 1.2 cP. The relative nondimensional elastic parameter is the elastic capillary 
number, Cac = /J,ja/G m where \i is the fluid viscosity, 7 is the local shear rate, a a length scale 
based on the particle size, and G m is the elastic shear modulus. 
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Hagen— Poiseuille of Cellular Blood Flow 

In these simulations, the capillary number of the particles increases with the radial distance from the 
axis of the tube, reaching a maximum at the tube wall. The pressure-gradient for these simulations 
is constant and is mimicked by a body force in the axial direction. Periodic boundary conditions 
are implemented at the tube inlet and outlet to allow the flow of RBCs exiting the domain to re- 
enter. The fluid domain for these simulations is 128x128x512 (8.39 million grid points). Clusters 
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of RBCs form in relatively high volume fraction Hagen-Poiseuille flow, and these clusters increase 
the relative suspension viscosity. Statistics captured include the number of clusters, average cluster 
size, and the effect of hematocrit and capillary number. The formation of three of these clusters 
is visualized at different time steps for the simulation described here. In Hagen-Poiseuille flow, 
the shear stress is at a maximum at the tube wall. As a result, the RBCs nearest the tube wall 
experience the highest rate of shear. This relatively high shear results in higher capillary numbers 
for the RBCs near the tube wall. The particles also tend to migrate away from areas of high shear 
and form a core in the center of the tube where the RBCs remain relatively undeformed. 

Linear Shear Flow of Red Blood Cells and Capsules 

Both wall-bounded and unbounded shear simulations are presented in this movie. The RBC simu- 
lation is in a cubical fluid domain of 256x256x256 (16.8 million grid points) with 2470 deformable 
RBCs (40% hematocrit) at Ca^ = 0.02. The movie portrays the different deformed shapes of the 
biconcave RBC geometries [5]. The unbounded shear simulation is of deformable spherical capsules 
(Cac = 0.02) with RBC membrane properties in a 144x144x144 domain. The video shows clus- 
ter formation along the compressional axis and cluster separation along the extensional axis with 
corresponding changes in stress on the particle's surface. 
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